Matter waves in the Talbot-Lau interferometry. 
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Paths of particles, emitted from distributed sources and passing out through shts of two gratings, 
Gq and d, up to detectors, have been computed in details by the path integral method. The slits 
are represented by Gaussian functions that simulate fuzzy edges of the slits. Waves of matter be 
computed by this method show perfect interference patterns both between the gratings and behind 
the second grating. Coherent and noncoherent the distributed particle sources reproducing the 
interference patterns are discussed in details. Paraxial approximation stems from the wave function 
when removing the distributed sources onto infinity. The more hard-edged slits of the grating Gi 
are examined by simulating those slits by a superposition of more hard the Gaussian functions. 

As for the particles here we consider fullerene molecules. De Broglie wavelength of the molecules 
is adopted equal to 5 pm. 

PACS numbers: 03.75.-b, 03.75.Dg, 42.25.Hz 



I. INTRODUCTION 

Interferometry with matter waves, where par- 
ticles are presented by heavy molecules, such 
as fullerene molecules, attracts the last years a 
great interest of scientific community; see review 
of Cronin, Schmiedmayer, and Pritchard [Ij, and 
rich set of references ibid. Interferometric exper- 
iments help to disclose the very basic principles 
of quantum physics with systems of rather large 
size and complexity [2HQ. 

Heavy molecules, having masses about 100 
amu and more, are particles showing under or- 
dinary circumstances almost classical behavior. 
Indeed, diameter of the fullerene molecule Ceo, 
see Fig. [l| is about 0.7 nm [9j, but de Broglie 
wavelength is ^ 5 pm [6l [lO] (velocity of the 
molecule is about 100 m/s). The wavelength is 




FIG. 1: The fullerene molecule CgQ consists of 60 carbon 
atoms. From the left, a characteristic size of the molecule to- 
gether with its characteristic de Broglie wavelength are shown. 

shorter of the diameter by a factor of 100 and 
much smaller than distances between atoms in a 
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solid body, whereas a size of the molecule exceeds 
these distances. Interference fringes emergent 
at scattering such heavy particles on a grating 
bear an important information about the wave- 
particle duality of these particles. Therefore, 
in order to reveal wave-particle duality for the 
heavy particles, interferometers should possess 
by heightened requirements |7j. 

By adopting the wave-particle duality formal- 
ism relative to the heavy particles, one may then 
deal with fiows of such particles as waves inci- 
dent on the slit grating [TTHllj. Methods based 
on computing the Fourier images of general two- 
dimensional periodic objects in light and elec- 
tron optics [inHIHj can be proposed as good algo- 
rithms for finding interference patterns induced 
by the matter waves. Next, a method adopted 
for the analytical description of diffraction uses 
the Kirchhoff-Fresnel integral [6]. The integral 
represents spherical waves emitted from a radi- 
ation source. Because the integral is taken only 
over the slit, the slit appears to act like an ef- 
fective source of radiation, in accordance with 
the Huygens principle [19] . This principle states 
that wave fronts passing through slits act as ef- 
fective sources, reradiating spherical waves. 

It is not intuitively obvious that the particles 
passing through the slits should be defiected to 
recreate spherical wave patterns. Paradox is that 
how does the particle identify positions of all 
slits? Nairz, Arndt, and Zeilinger [8j have writ- 
ten "The wave-particle duality states that the 
description of one and the same physical object 



suggests the local particle picture in the source 
and on the screen, but a wave model for the un- 
observed propagation of the object." While the 
particles go through the slits one by one, a single 
particle always gives a single click at a detector. 
Between the source and the detector history of 
the particle remains unknown. We may consider 
all paths of point-like particles going through the 
slits and converging into the detector. This task 
is solved by the path integral method [20.^ ^Ij. 
The converging paths at the detector give rise to 
superposition of different histories, which pro- 
duce an interference phenomenon. 




FIG. 2: The Talbot-Lau interferometer scheme: two gratings, 
Go and Gi are situated in consecutive order along a particle 
beam emitted from a distributed source. 



Here we will consider emergence of the interfer- 
ence patterns at passing fuUerene-like particles 
through two gratings, Go and Gi, placed in con- 
secutive order along a particle beam, see Fig. [2} 
It is a typical scheme of the Talbot-Lau inter- 
ferometer [^ [23] . For finding the interference 
pattern we will compute all paths going from 
the source to the detector screen, i.e., we will 
compute the path integrals. Emergence of the 
interference patterns is considered in the near- 
field - distance between the grating is half of the 
Talbot length, Zt/2 = d'^/XdB = ^i — ^o- Here 
AdB is de Broglie wavelength and rf is a distance 
between the slits. 

Coherence properties of the beam are indis- 
pensable conditions for manifestation of the 
wave-particle duality and, as an effect, emer- 
gence of interference fringes. These properties 
depend on the source, collimation slits, and other 
extra conditions [8j . Ones distinguish spatial co- 
herence and spectral coherence. The first is con- 
ditioned by finite width of the source that can 
exceed de Broglie wavelength. And the second 



can deteriorate because of generating the parti- 
cles with different escaping velocities. Here we 
will try to reproduce the both coherence proper- 
ties. We will consider partially coherent sources, 
referred to as the generalized Gaussian Schell- 
model (GSM) sources [241427] . Such sources re- 
produce planar GSM beams. These beams are 
characterized by two signs, that are an inten- 
sity distribution across the beams and a spectral 
degree of the coherence. The both are repre- 
sented by the Gaussian distributions with their 
own dispersion constants aj and ag. With suit- 
able choices of aj and a^, such a source will sim- 
ulate a beam, called GSM beam [24> .25]. 

Observe, however, that the particle beam rep- 
resents, in accordance with the particle- wave du- 
ality, a matter wave, that spreads through the 
grating device. It should be noted, that the wave 
consists of two parts, amplitude and phase ones 
multiplied each other, |^) = y/p exp{iS/h}. 
Only on the detector we register the intensity 
proportional to p = (^|^). We will apply the 
Gaussian Schell-model to the wave function as a 
simple Gaussian form-factor loaded by either of 
two dispersion constants a/ and ag. 

The article is organized as follows. In Sec. H 
we compute passing the particles through the 
two gratings, Go and Gi, being emitted from 
distributed sources. Exact formulas describing 
interference effect emergent both between the 
gratings and behind the second grating are rep- 



resented in this section. Sec. [Hi] represents cal- 
culations of interference patterns emergent at 
scattering fuUerene-like particles on the gratings 
with finite amount of slits. Here we consider an 
effect of spatial coherence to interference pat- 
terns. In Sec HVl we remove the sources to in- 
finity and obtain the formulas for interference in 
the paraxial approximation. Here we consider 



an effect of spectral coherence. Sec. \V\ discusses 
interference from the gratings with more hard- 



edged slits. Sec. VI is concluding. 



II. 



PATH INTEGRAL: PARTICLES PASSING 
THROUGH SLITS IN TWO SCREENS 



Computation of passing a particle through the 
system of two gratings. Fig. [2} is based on the 



path integral technique [211 EH]- Let us begin 
with writing the path integral, that describes 
passing the particle through slits made in two 
screens. They have been preliminarily prepared 
in two opaque screens situated perpendicularly 
to axis z, see Fig. (S) For this reason we need to 
describe a movement of the particle between the 
screens and possible changing in its deflection at 
crossing the slits. We believe, that between the 
screens the particle moves as a free particle. Its 
Lagrangian is 



m — - + const. 



(1) 



Here m is mass of the particle and x is its 
transversal velocity. A longitudinal momentum 
Pz is much greater than its transverse compo- 
nent [HI [m [I3j and we believe it is constant. 
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FIG. 3: Passage of a particle along path Xs ^ xq ^ xi ^ X2 
through two screens containing slits with widths 26o and 26i . 
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The particles flying through the flrst slit along 
a ray a will pass this slit, ranging from xq — 6o 
to Xo + 605 almost surely. And let its deflection 
on the flrst slit be such, that it passes through 
the second slit ranging from Xi — hi to Xi + 61. 
The path integral in that case reads 



?/;(x2,Xi,Xo,Xs) 

hi 

j K{x2, r + To + ri; Xi + ^1, T + tq) 

-hi 
ho 

j K{xi + ^1, T + To; xo + ^0, T) 



-ho 



xK{xq + Co, T; Xs, 0)d^od^i. 



(2) 



Integral kernels (propagators) for the particle 
freely flying are as follows [21j 
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By substituting the propagators into the path 
integral (J2]) we get 
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The both integrals we can see are computed 
within flnite intervals [— 605+&0] and [—61, +61], 
respectively. Observe, that the integrating can 
be broadened from — oc to +oc. But in this 
case we need to load the integrals by the step 
functions equal to unit within the flnite inter- 



vals [— 605+&0] and [—61, +61] and they vanish 
outside of the intervals. 

In order to obtain exact solutions of the path 
integral let us take supposition, that the slits 
possess slightly fuzzy edges. In that case, the 
form factors, simulating the slits, can be repre- 



sented by the Gaussian function \21\ 

GiO = exp{-e/2b'}. (5) 

Effective width of the Gaussian curve is given 
by the parameter b. About two thirds of area of 
the curve are placed between points —b and +6, 
and one third of the area is beyond the interval 
(—6, +6). Consequently, tunneling through the 
screen in the vicinity of the slits we believe can 



take place. It is due to the fuzzy edges adopted 
above. 

Let us insert into the integrals the Gaussian 
factors G{^o) and G(Ci) with the parameters bo 
and bi and replace the finite limits ±6o and ±6i 
by the limits from — oc to +oc. Take into ac- 
count, that we have a solution for the inner inte- 
gral, that has been written out in [29j. Now we 
can rewrite the integral Q more definitely 
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>d^i. (6) 



Solutions of such integrals stem from the formula In order to get such a form, let us regroup all 

terms under the integral ([g]) and collect their at 



/ 
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-l3^/4a^-f 
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(7) 



coefficients ^f, ^i, and free from it. We have 



1. the term at ^f: 



im ( 1 1 ih 
Oi = -r\— + — + 



2h\Ti To mhl T'^{l/To + l/T + m/mhl) 



(8) 



2. the term at ^i: 



/5 = -2 



im ( {x2 — xi) (xi — xo) 



2h 



Tl 



To 
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il/ro + llT + m/mhl) j' 



(9) 



3. the term free from ^i: 



im I {x2 — xiY (xi — xoY 

2h \ Tl To 



{xo - XsY {{xi - xo)/to - {xo - Xs)/Tf 



T 



[1/to + 1/r + ih/ml 



(10) 



Now we can express the terms (7r/(— a))-^/^ and 
7 — /3/4a in the right part of Eq. ([7|. The first 
term, accurate to the multiplicand 
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relates to an amplitude factor A of the wave 
function '0(^2, Xi, Xq, Xg). For this reason, by 
multiplying the term (7r/(— a))-^/^ by the mul- 
tiplicand (11) we obtain the amplitude factor A. 



It has the following view 
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2TrihT 
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ihTi 
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In turn, the term 7 — /3^/4q; reads: 
7 - /3^/4a 



2h 
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mbl{l + To/T)J To 



{X2 - Xif (xi - Xpf {xo - Xsf {{xi - xo)/to - (xq - Xs)/Tf^ 
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To 
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(12) 



(13) 



A. Series of replacements 



Let us now define effective slit's half- widths 



^.7,0 






(14) 



where j takes numbers and 1. And next we 
determine a complex time-dependent spreading 



G 



j^^j 



^j,o + i 



hTj 



2mcFn 



(15) 



where j = 0, 1 and at j = we suppose r_i = T. 
It should be noted, the complex spreading pa- 
rameters have been presented in works [30^ ^J . 
Here we choose the same representation. 

More one step is to replace fiight times T, tq, 
and Ti by fiight distances (zq — ^s), (^i — ^o), and 
(^2 — ^i). Fig. |3J This replacement reads 

T = (zo - z^)/v^, 

To = {zi - zo)M, (16) 

^ Tl = {Z2 - Zi)/V^, 

where Vz is a particle velocity along the axis z. 
We note that mVz = Pz is z-component of the 
particle momentum. Next, we introduce the de 
Broglie wavelength A^^ = h/pz^ where h = 2tiTi 
is the Planck constant. Rewrite the complex 



spreading (15) according to these remarks 
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(17) 
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Hereinafter, for brevity, we will not write the 
subscript dB at A. 

Define a dimensionless complex distance- 
dependent spreading as follows 



^J>j - 
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^j-i 
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^j) 



4^^j\o 



Also we define a dimensionless parameter 

(xo - Xs) {zi - Zo) 



--0 



1 



(18) 



(19) 



(^0 - ^s) (^1 - ^o) 

that tends to 1 as Zg ^ — oc. Now we can 
rewrite the terms A and 7 — /3^/4Qf, represented 



in Eqs. (12) and (13), by replacing cumbersome 



parameters by defined above. 
The amplitude factor, ^4, of the wave function 

'0(x2, Xi, Xo, Xg) has the following view 

A=,[Z^. ^,^ ^„ . . (20) 



27rifir i^(So,,o,Si,,J' 

Here we do not replace T by (zq — z^jvz. It 
points out to distance to the source accurate to 
division by Vz^ As we remove the sources to infin- 
ity (T -^ oc) the amplitude factor tends to zero. 
Nevertheless, we keep a finite value of the pa- 
rameter A as a normalization factor of the wave 
function [29j. 

The phase term, 7 — /3^/4a for the same wave 
function '0(x2, Xi, Xq, Xg) reads 



7 - /5^/4q; = ill 



' {X2 -Xif {Xi - XqY 

A(^2-^i) Mzi-zo) 
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r 



(21) 



The term D(So,zo5 ^1,^1) represented in divisors 
of expressions ([20| and ([21]) is as follows 



(22) 



Thus, we will have in mind that the expres- 
sions (21) and (22) contain full information 



^(So,2:o5 Sl,2;i) — \^0,zo^l,zi " 



^2 - Zi 
Z\ — Zq 



B. Matter waves behind the gratings Go and Gi 

For observation of the wave field behind the 
second grating, Gi, we need to situate a detec- 
tor in a point (^2, ^2), Fig.[3J On the other hand, 
if we wish to observe the wave field between the 
first and second gratings, between Gq and Gi, we 
must situate the detector in a point (xi,Zi). In 
order to realize the second case, it is sufficient to 
put in the expressions (21) and (22) X2 = Xi and 
Z2 = Zi. Let now the point (^2, ^2) = (^i, ^1) be 
situated in a region between the gratings. Ob- 
serve, in this case, that all terms containing dif- 
ferences (^2 — Zi) vanish. Next, we move the 
detector in the point (xi,Zi) -^ (x^z) that can 
be situated anywhere between the gratings. 



about wave fields both between the gratings and 
behind the second grating. It depends on choos- 
ing coordinates of position of the detector, either 
Z2 -^ z^ X2 ^ X, or Z2 = Zi ^ z, X2 = Xi ^ X. 
We write variables (x, z) instead of (^2, ^2) at de- 
scribing the wave field behind the second grating. 
And we write the same variables (x, z) instead of 
(xi,zi) when we describe the wave pattern be- 
tween the gratings. In the last case, we have 
replacements (^2 — ^1) -^ {z — Zi) ^ (z — z) = 0, 
(x2 — Xi) -^ {x — Xi) ^ (x — x) = 0. Terms con- 
taining such differences either vanish, or become 
unit if a ratio {x — x)/{z — z) takes place. Now we 
can write out the wave patterns at passing the 
particle through single slit within each grating. 
A particle emitted from the point {x^^z^)^ be- 
longing to the source, can be detected in the 
point (x2,^2) ^ (^5^)5 z > Zi (behind the sec- 
ond grating) in accordance with the following 
wave function 



ip\X^ Z, Xi, Xq, Xg) 
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27vihT 
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(23) 



In the region between the gratings, z < Zi, we employ replacement (xi, Zi) ^ (x, z) and consequently 



?/;(x,z,Xo,Xs) = 



m 



^^^i^^^^-^^T" 



{X - XqY [^ _ ^|_^ ^ (xq - XsY 



X{z - zo) \ Eo, 



^0 , 
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(24) 



As we stated above, the function '0(x, z, Xq, Xg) stems from '0(x, z, Xi, Xo,Xs) as soon as we put 



in Eq. ([23]) values Z\ = z and X\ = x. Observe, 

that L^(So,^o,Si^^J = T}^'^ ^ so far as Si^^, = 1 at 
z = Zi = ^2, see Eqs. (18) and (22). 



III. INTERFERENCE PATTERNS FROM 
TWO-GRATING STRUCTURE 

The gratings shown in Fig. [2] we beheve have 
different amount of shts. The first grating has 
A^o shts, so no = 0, 1, • • • , (A^o — !)• The second 
grating has A^i shts, and ni = 0, 1, • • • , (A^i — 1). 
Distance between shts for the first grating is rfo 
and for the second grating is di. Now we can 
write out a wave field beyond the gratings that 
is composed of superposition of wave functions 
representing coherent emission from all slits. 

The matter wave in a zone between the first 
and second gratings reads 






(25) 



And the matter wave extending beyond the sec- 
ond grating reads 



|^i(x,z,Xs,A)) 



J—1 _~ I f 

Y^ ^ i)\x,z, Ini - 



A^i-l 



ni=0 no=0 



rfl, 



no 



A^o-1 



rfo,Xs . (26) 



The both wave functions are represented without 
normalized factors. We need no here in these 
factors, since our interest is to show a general 
pattern of the density distribution 

p(x,z) = (^(x,z,Xs,A)|^(x,z,Xs,A)). (27) 

As for the particles supporting the matter 
wave here, we consider fuller ene molecules [lOj. 
They are massive molecules. Its mass is about 
TTiCeo ^ 1-2 X 10~^^ kg. The fullerene molecule 
has radius i^Ceo = 350 pm [9j. In turn, the 
de Broglie wavelength is much smaller than 
this radius. For instance, the de Broglie wave- 
length equal to 5 pm appears in the experimen- 
tal work [10]. The wavelength AdB = 5 pm 



is adopted in this work as well. The widths 
of the slits (open slit windows) are as small as 
26o = 75 nm in Go and 2h\ = 150 nm in Gi. 
Distances between the slits are equivalent for the 
both gratings do = di = 500 nm. The distances 
have been double increased in contrast to those 
given in [lOj and [29]. It is done in order that 
we could see definitely the Talbot carpets aris- 
ing between the gratings, and maybe beyond the 
both gratings. 

Here we consider example of emergence of in- 
terference patterns from the gratings Go and Gi 
containing even amount of slits, A'o = 32, and 
odd amount of slits, A^i = 33. Difference in par- 
ity of the numbers A'o and A^i is conditioned by 
the fact, that the first self-image of Go is shifted 
exactly on half period of the grating. Due to 
this trick, the slits of Gi are located exactly on 
nodes of the first self-image of the grating Gq. 
The grating Gi keeps "open gates" for particles 
to spread further. 

Fig. [4] shows the density distribution p{x^ z) in 
the near-field region for different positions of a 
point source: (a) Xg = /xm; (b) Xg = 2 /xm; 
and (c) Xg = 4 /xm, and at z^ distant from the 
first grating on —0.5 m. Here we have shown 
cases of shifting the point source to a positive 
area, Xg > 0. As for negative shifting it can be 
obtained by simple refiection of the interference 
patterns about the plane (^, z) intersecting axis 
X ai X = 0, see frame of axis in Fig. [2} That 
is, the slits should be subjected to the following 
inversion 12 ^ 19, 13 ^ 18, 14 ^ 17, 15 ^ 16. 

We can see in Figs.[4](a), [4](b), and|4](c) ripples 
of high order against the background of the basic 
divergent rays. They are induced by presence 
of many lateral slits invisible in these figures. 
The interference patterns are seen to change at 
changing position of the point source. 



A. Coherent and noncoherent sources 

Let us now turn to the coherence properties of 
the particle beam [8j. We will assume, that point 
sources are coherent, if wave functions relating to 
these sources are summed together in the vicin- 
ity of a detector. And they are noncoherent if in- 
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FIG. 4: Density distribution pattern p(x, z) in the near- field 
region z G (0, 1.5zt) = (0,0.15) m, ATq = 32, iVi = 33, de 
Broglie wavelength Ads = 5 pm, Talbot length zt = 0.1 m. 
Distance from Go to the source is Zs = —0.5 m : (a) the source 
is situated on the optical axis, Xs = 0; (b) the source is shifted 
from the central axis on Xs = 2 /xm; (c) the source is shifted 
from the central axis on Xg = 4 /im. 
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FIG. 5: Density distribution pattern p{x, z) averaged over 
all point sources localized with increment 6xs = 0.25 /xm in 
the interval Xs G [—4,4] /xm, Zs = —0.5 m, A^o = 32, Ni — 
33. The averaging has been carried out with the Gaussian 
kernel (29) loaded by the dispersion constant ai\ (a) coherent 
beam, Gj — 10 /xm; (b) almost coherent beam, a^ — 1 /jui; 
(c) almost noncoherent beam, a^ =0.3 /xm. 



tensities (probability densities) are summed. In 
particular, superposition of the wave functions in 



the vicinity of the detector can be with weakened 
contribution of cross terms, that provide inter- 



ference effects. This case relates to intermediate 
variants, which can be considered by appeahng 
to the Gaussian ScheU-model [SUES]. We have 
to keep in mind, however, that the wave func- 
tions are primary quantum subjects, whereas in- 
tensities are found on a final stage at reading 
from detectors. 

Averaging of the wave functions from all point 
sources distributed along Xg at fixed Zg is carried 
out in the following form: 

P{X,Z) = (Tj^Yl 

I II 

Xg Xg 

{^{x,z,x[,\)\^{x[,xl ,a ^)\^{x,z,xl ,\)) . (28) 
The Gaussian kernel ^{x^,x^ , gj) reads 

1 r {x[-xy- 



^[Xg, Xg , CTj-J 



2na, 



exp 



2 a/ 



(29) 



with the dispersion parameter a^ being an ef- 
fective coherent width of the beam. This pa- 
rameter presented as a factor in front of the 



sums (28) provides identity of dimensionalities 



for the probability density distributions. 

First, one can notice that at a^ <C 1, the Gaus- 
sian kernel fits the Dirac (^-function. And the 



expression (28) drops to a simple summation of 



the probability densities 

p{x, z) - ^(^(x, z, x,, A)|^(x, z, x,. A)). (30) 

Xg 

And at cr^ :» 1 the Gaussian kernel degenerates 
to a constant. In that case we have 

p{x, z) - ^(^(x, z, x,. A) |^(x, z, X3, A)) 

Xg 

+ E E(*(^'^'^s'^)l*(^'^'^s",A)). (31) 

I 1 II II 

The second sum here contains the cross terms, 
that introduce interference effects from different 
point sources. 



The expression (30) represents an example of 



wholly noncoherent beam. Whereas, the expres- 



sion (31) gives completely coherent beam. Inter- 



mediate coherence beams are possible as well. 

Figs. [5](a), [5](b), and[5](c) show averaged den- 
sity distributions for different depth of coherence 
ranging from a coherent beam to noncoherent 



one. Summation is taken for all point sources 
localized at Xg = —4, —3.25, —3.5, • • • , 4 /xm, the 
increment is 5xs = 0.25 /xm. Distance to the 
source is Zg = — 0.5 m. The figures demonstrate 
interference patterns that are reproduced from 
(a) coherent beam, g^ = 10 /xm, to (c) almost 
noncoherent beam, a^ = 0.3 /xm. The high-order 
interference fringes are seen to be washed out 
as the dispersion parameter a^ decreases from 
10 /im to 0.1 /xm, what is in good agreement 
with computational results given in 

5' 

(a) 




x.jum 



0.8 

H 

0.4 




(b) 



-1 



x,//m 



1 



FIG. 6: Interference fringes observed in the cross-section of 
the density distribution pattern (28) at z = zt'- (a) almost 
coherent beam, cr^ = 1 /xm; (b) noncoherent beam, a^ = 
0.1 /im. Gray vertical strips indicate opaque spaces between 
slits in the grating Gi . Arrow points out a level of pedestal. 

Let US put a detector screen on the distance 
L = Zt from the grating Gq and look out 
on emergent interference fringes for two cases, 
namely, for a^ = 1 /xm and for a^ = 0.1 /xm, 
see Figs.|6](a) and[6](b). One can see, that at de- 
creasing the dispersion parameter a^ a pedestal 
supporting the interference fringes emerges. Po- 
sition of the pedestal is pointed out by arrow in 
Fig. [6) The pedestal can be found as an absolute 
minimum for all the interference fringes 

Pmin = minp(x,z^) (32) 

We can find also an absolute maximum 

i^max = max p(x, z^) (33) 
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The interferometric visibility be computed by a 
formula 



V 



+ Pu 



(34) 



quantifies contrast of the interference fringes. 
_ (a) pedestal 
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(b) visibility 
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FIG. 7: Pedestal Pmin (a) and visibility V (b) calculated as 
functions of the effective coherent width a^ of the beam. 

Variations of two parameters, level of the 
pedestal Pmin and the visibility V^ simulated as 
functions of the effective coherent width of the 
beam a^ ranging from 10~^ /xm to 10^ /xm are 
shown in Figs. [7|(a) and [7|(b). A crossover is 
clearly seen within the interval 1 /xm down to 
0.1 /xm. It represents a smooth transition from 
the coherent source to noncoherent as the dis- 
persion parameter a^ decreases within the men- 
tioned interval. 

Surprisingly, if radiation from all point 
sources, all situated near the interferometer, is 
coherent, superposition of these radiations repro- 
duces a perfect interference pattern as if from 
a single remote source, see Fig. [5](a). At that, 
all point sources situated near the grating Gq 
demonstrate radial divergence of rays beyond the 
grating Gi, as is seen in Figs. |4](a)-(c). To see 
a general pattern of such a radial divergence we 
have simulated emergence of an interference pat- 



tern from two gratings consisting of A^o = 32 
and Ai = 33 slits. Fig. [8](a) shows this inter- 
ference pattern. Because of closely spaced the 
single source, z^ = —0.5 m, we see, that there 
are three groups of divergent rays from the sec- 
ond grating. Lateral rays of these groups inter- 
fere on a distance of the half Talbot length from 
the second grating. It is remarkable, that su- 
perposition of these divergent rays from differ- 
ent coherent sources reproduce the Talbot-like 
interference pattern as shown in Fig. [5](a). 

Fig. [8](b) demonstrates the very interference 
but for the point source on Zg = — 50 m distant 
from the grating Gq. One can believe in this 
case, that the source has been removed almost 
to infinity. One can see, that Fig. |8](b) displays 
in the near-field the Talbot carpets, that con- 
tinue behind the second grating. For emergence 
of the perfect Talbot carpet it is necessary to 
satisfy the following three requirements J33H35] : 
(a) a particle beam is paraxial; (b) ratio of de 
Broglie wavelength, AdB, to period of a grating, 
rf, tends to zero; (c) number of the slits tends to 
infinity. Fig.[9]shows the Talbot carpet emergent 
from the two gratings configuration consisting of 
Ao = 64 and Ai = 63 slits. The Talbot carpet 
looks the better, the more number of the slits is 
in the gratings, ideally tending to infinity [36j. 

It is instructive to compare patterns shown in 
Figs. [5](a) and |9| The first Talbot-like pattern 
was got by superposing many coherent rays aris- 
ing from sources situated near the interferome- 
ter. This pattern is washed out as soon as the 
rays become noncoherent. Whereas the second 
Talbot pattern is got from a single plane wave 
incident to the interferometer from infinity. 

Observe that in case of the coherent monochro- 
matic beam, presence of the second grating does 
not affect on emergent the interference pattern. 
That is, the same interference pattern emerges if 
we would remove the second grating. 

Briefiy, the incident beam shows spherical 
equiphase surfaces of the matter wave when the 
source is positioned nearby the interferometer. 
At removing the source onto infinity, the spher- 
ical equiphase surfaces degenerate into planar 
equiphase surfaces. In that case the spherical 
wave turns into the plane wave incident on the 
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FIG. 8: The density distribution pattern p{x,z): (a) particle source is spaced closely to the first grating, Zs = —0.5 m; 
(b) particle source is spaced far from the first grating, Zs = —50 m. De Broglie wave length AdB = 5 pm, and zt = 0.1 m. 



^dB = 5 pm 




FIG. 9: Talbot carpet on two grating configuration: the 
gratings Go and Gi consists of A^o = 64 and Ni = 63 slits, 
respectively. De Broglie wavelength Aas = 5 pm, distance 
between slits d = 500 nm, the Talbot length z^ =0.1 m. 



pendicular to the equiphase surface and particles 
pass in parallel to the axis of the optical system, 
the axis z. Such a particle beam is called the 
paraxial beam. In the next section we obtain 
this paraxial approximation. 

IV. PARAXIAL APPROXIMATION 



The wave functions (23)- (24) will describe in- 



terference in the paraxial approximation as soon 
as a limit z^ ^ —oo will be reached. Observe, 
that in this limit Sq = 1, ^o,zo = (^o,zo/(^o,o and 



1- 



1 



So. 



zo 



A(^l - Zq) 
47rcro,0 CTo.zo 



(35) 



interferometer. The particles momenta are per- For verification see Eqs. (17)-(19). 
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Next, we need also to reinterpret the term 



^(c^O,zo,Cri^^J 



z — Zq \ ai^zi ^0, 



ZQ 



Z — Z\ 



I. (36) 



\ \^\ - ^0/ CTi^o Cro,o V^l ~ ^0, 

It is instructive to compare this expression with 



Eq. (22) representing the term i^(So,zo5 ^1,21)- 

Here and in the next formula signs colored in 
blue relate to an area reaching out after the sec- 
ond grating. Remaining signs colored in black 
deal with an area between the gratings. 

As soon as the all reductions have been done, 
we obtain the wave function in the paraxial ap- 
proximation, i.e., with the source remote onto 
infinity. It has the following view 



'0(x,Z,Xi,Xo) 



A., 



-jexp{i.[( 



' [X - Xi)^ . (Xi - Xq) 

\{z-zx) 



D{ao^zQ,ai^ 
X{z - Zi) CTq,^, ^ 



47rao,ocro,zo 
2 



(x - Xi) . (xi - Xo) 



^0. 



(37) 



Here a factor A^o replaces Jm/{27rihT). Since 
at T ^ oc the quadratic root tends to zero, then 
^4^0 tends to zero as well. We will ignore this 
fact, and suppose let ^4^0 be some constant. 



The wave function (37) describes interference 



effects emergent both behind the second grating 
Gi and between the gratings Go and Gi . In order 
to get the wave pattern between these gratings it 
is sufficient to set Zi = z and Xi = x. In this case, 
ai^zi becomes ai^o and Djap^zn^ ^i.z^) reduces to 
(c'"o,zo/^o,o)^^^- All terms in (37) containing dif- 



ferences (x — Xi) and {z — Zi) disappear. To put 
it bluntly, all terms in (36) and (37) colored in 
blue should be removed. 

Let us look on interference patterns emergent 
on such a device at illumination by matter waves 
with different wavelengths. Figs. [Tol(a) to [Tol(c) 
show interference patterns for cases of incident 
particles having different de Broglie wavelengths: 
AdB = 3 pm, AdB = 5 pm, and AdB = 7 pm. 
Respectively, velocities of the fuUerene molecules 
at given wavelengths are v^ ?^ 184 m/s, v^ ^ 
110 m/s, and v^ ^ 79 m/s. 

The grating Gi is situated at half of the Tal- 
bot length. Given AdB = 5 pm and d = 500 nm 
the Talbot length is Zt = 2(i^/AdB = 0.1 m. 
Depending on the de Broglie wavelength chosen 
the interference pattern, emergent between the 



gratings, discloses different scaling. Because of it 
different interference patterns behind the second 
grating are formed. Most intensive the interfer- 
ence pattern arises at the de Broglie wavelength 
equal to 5 pm, since position of the grating Gi 
has been tuned on the first self-image of Gq aris- 
ing at the same wavelength. 

One can observe a resonance effect at crossing 
the first self-image of Go by the grating Gi. It 
can be achieved by changing the wavelength A^^ 
at crossing the resonance condition A^^^ = 2rf/z^ . 
In the case under consideration A^^^ = 5 pm. Ob- 

dB ^ 

serve, that a maximal emittance from the grat- 
ing Gi is at A^j3 = A^^^, when it is positioned 
exactly on the first self-image of Go- There can 
be also high harmonics at nX^^ = A^^^ (here n is 
integer), when Gi is positioned on the high or- 
der images of Gq. In these cases the emittance 
from Gi quickly drops off with increasing n. The 



emittance to be expressed by a parameter (33) is 



shown in Fig. [11} Instead of representing via de- 
pendence of the de Broglie wavelength, here we 
show the dependence via the fuUerene velocity 
v^ = h/{m^ \b)- Here h is the Planck con- 
stant. Evaluations say that for velocities 100- 
250 m/s one needs to support ultralow temper- 
atures ranging about 3 • 10~^ to 3 • 10~^ K. 

Let us suppose that the remote source emits 
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FIG. 10: Density distribution pattern p{x, z) at different 
de Broglie wavelengths: (a) A^^ = 3 pm, v^ ^ 184 m/s; 
(b) A^B = 5 pm, v^^^ ^ 110 m/s; (b) A^^ = 7 pm, v^^^ ^ 
79 m/s; No = 8, Ni = 9 and Talbot length zt = 0.1 m. 



particles with different de Broglie wavelengths. 
Distribution over all wavelengths submits to the 
Gaussian with average AdB = 5 pm and disper- 
sion constant ag = 2.25 pm. An averaged inter- 




T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1— 

100 200 Vp , m/s 300 

FIG. 11: Emittance Pmax in the cross-section z = z^ vs 
velocity of the fullerene molecules. Arrows 1 and 2 point out 
to the first and the second resonance harmonics. 



H 



FIG. 12: Density distribution pattern p{x, z) at A^o = 8, 
Ni — 9 from a matter wave containing particles with different 
wavelengths from A^g = 2 pm to A^g = 8 pm. 



ference pattern for the wavelengths ranged from 
3 pm to 8 pm with increment 5\^b = 0.25 pm, 
and under assumption that the sources are non- 



coherent, is shown in Fig. 12, One can see, that 
fine-structured details in the interference pat- 
tern disappear. They are simply washed out. 
Nairz, Arndt, and Zeilinger [8j have written for 
that occasion: "because the detector records the 
sum of the correspondingly stretched or com- 
pressed diffraction pictures, the interference pat- 
tern would be washed out. And in contrast to 
the spatial contribution, there is no gain in longi- 
tudinal (spectral) coherence during free fiight." 



Fig. [121 confirms aforesaid thought. Interference 
patterns disclose equivalent image blur irrespec- 
tive of choosing of the dispersion constant ag. 
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XQ+b X 



FIG. 13: Approximation of the step function by set of the 
Gaussian functions presented in (38) with (a) 77 = 1, K = 8; 
and (b) r] = 1.5, K = 16. 



V. GRATINGS WITH MORE HARD-EDGED 
SLITS 

Let us compute the path integral ^ for case 
with more hard-edged shts. With that aim we 
should fill the slits uniformly by a number of 
the Gaussian functions ^ with more sharp bell 
curves. The step function, that simulate a single 



slit, can be approximated by the following a set 
of the Gaussian functions 



Gi^,b,rj,K) = 




K 



^exp< 



k=i 



{Kj - h{K - {2k - i))y 
2(6?7)2 



.(38) 



Here parameter 6 is a half- width of the slit, real 
77 > is a tuning parameter, and K can take 
integer values. At K ^ oc this function tends 
to an infinite collection of the Kronecker deltas, 
that fill everywhere densely the step function. 
Fig. [T3I shows approximation of the step function 



by the set of the Gaussian functions (38) with 



(a) 77 = 1, X = 8 and (b) r] = 1.5, K = 16. 

The form- factor G(^i) in Eq. ^ is replaced 
here by the function G{^i^bi^r]i^ Ki). Solution 
of the integral ([6]), containing the set of the 
Gaussian functions (ISSl), satisfies the formula 




,c^fcC^+/5fcC+7fe 



d^ 



^-/3|/4afc+7fe 



(39) 



Before we compute parameters ak^ /3k^ Jk '^^ ^^ 



necessary to open the square in Eq. (38). After 
all we find these parameters 



1. the term at ^f: 



2. the term at ^i: 




2\ / ' 



Pk = -2 



+ i 



2h \ti To mr]l bj tHI/tq + 1/T + ih/mbl) 



imf{x2-xi) (xi-xo) , {xi-xo)/tq - {xo-Xs)/toT 




To 



{l/To + l/T + ih/mi 



mrji bi 



{K^-{2k-l)) ; 



3. the term free from ^i: 

im f {x2 — XiY (a;i — Xo)^ {xq 



Ik = 



2h \ Ti 
h 



To 



+ 



Xsf {{xx - Xo)/to - {xo - Xs)/Tf 



T 



(1/to + 1/T + m/mbl) 



+ \^^{K,-{2k-l)Y\. 
mrji I 



(40) 



(41) 



(42) 
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It is instructive to compare these parameters 
with those presented in Eqs. ([s]), ([9]), and (10). 
Differences easily strike the eye. The amphtude 



A minor difference of the amphtude factors (12) 



factor (12) for our new task reads 



and ( |43| ) is conditioned by an additional factor 
(2/7r)~/7yi presented in Eq. (39). In particular, 
at r^i ^ 1.5 and Ki = 1 the amplitude factor (43) 



A 



1 



will be twice as little of the factor (12). The de 



m 



r]i^ 71 yi27rhT D{h^M) 



(43) 



nominator D(6o, hi) here has the following view 



D(KM) 



\ 



tq \ mr]fbi{l + Ti/To) 



,^t\ ^mbi{l + To/T) To- 



rn 



In turn, the therm 7^ — /5|/4«fc reads 

im 



7fc-/5t/4«fc = 



2n 



{X2 - XiY {Xi -XqY {xq - Xsf {{Xi - Xo)/to - {xq - X^) /Tf 



Tl 



To 



T 



{{To + T)/TTo + m/mhl] 



mrii 



'(a;2-xi) (xi-xo) , {xi-xo)/to-{xo-Xs)/T h Ki ' 

H . ^x ,r^ —n—, — 7^77 + 1 ^^— (-"■! - {2k - ijj 



Tl 



To 



ro((ro + T)/Tto + ih/mbi) mill hi 



Vo + T^ m 



1 



Ttq mhl) ToTiD{bo,bi) 



(45) 



A. Series of replacements 



El 2, in this task is rewritten as foUows 



In order to execute series of the replacements 



as in Subsec. p A| it is proposed to introduce the 
foUowing scahng change for the half-width 61 



riib 



Ki 



1^6', 



(46) 



And consequently, we will have in mind that 
the parameters ai^o, ^i,Tn ^i,zi defined in 
Eqs. (14), (15), 18) can contain the very terms 
Ki and rji. After all, the parameters So,zo ^^d 



^0,2:0 



^1 ~ ^s . A(zi — Zq) 



Si, 



Z2 



Zs 

zo 



z\ 



Z\ — Zq 



+ 1 



4vrao2o 



^T^<^lo 



(47) 
(48) 



Here ao,o and ai^o have the original forms pre- 
sented in (14). 



The phase term 7/^ — /5|/4a/e, that has a rather 
complex form (45), reads 



ITT 



(x2 - xif (xi - x^y f 



7fc-/5|/4afc 

\[Z2 — ZijSp^^^ 

D(So,^o, Si,^J2 y\[z2 - Zi) \{zi - Zq) 



A(^2 - Zi) \{zi - Zo) 
(X2 - Xi) {Xi - Xo) ( ^ 



^2 
^0 



So, 



^ {xo-x,f ^ .^{K^-(2k-l)f 



^0 , 



^0 



+ i 



-'0,zo , 



X{zo - Zs) 

Kr{Kr-{2k-l)y 
27rr]jbi 



27ir)l 



(49) 



Here i^(So,zo5 ^1,^1) is represented in (22), but with Si^^-^ to be loaded from (48). 
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Now we can write the wave functions describing appearance of a particle both between the gratings 
and behind the second grating. 

B. Matter waves behind the gratings Gq and Gi 



Evolving a particle from a single slit in Gq to Gi is described by the very wave function (24): 

(x - x^f 



ip[x^ Z, Xo, Xg 



/ m 
i27rfirS. 



0,zo 



exp< ITT 



\{z - Zo) 



1- 



(xq - x^y 



+ 



(50) 



Observe that equivalence of the wave functions (50) and 24) is due to the fact that we use the same 



grating Gq in the both cases. But behind the grating Gi, having more hard-edged slits, movement 
of the particle from a single slit in Gi is described by a wave function 



exp<! 



ilj{x,Z,Xi,Xo,Xs,Ki) = 
{x — XiY {Xi — XqY ( 



m 



i27rhT 



771 V TT L^(So,^o,Si,^J 



E 
fc=i 



\X{Z-Zi) X{Zi-Zq)\ So,^o/ Hzq-Zs) 



El \ , (xo-x,f ^ .(i^,_(2fc-l))^ 



2Trrii 



X{z - zi)Eo, 



ZO 



[X — Xi) 



[.Xl - Xq) 



So \ , .K^{K^-{2k-l)y 



D(Eo,2o,Si,^J2 \^A(z-zi) X{zi-zq) 



^0,zo , 



+ i 



27rrilbi 



(51) 



One can see, that with Ki = 1 the wave func- 
tion (51), accurate to the factor {2/tiY^^ /r]i^ is 
simply the function (23). 



Interference patterns arise as superpositions of 
the above wave functions radiated from all slits 
of the grating Gq and from all slits of the grating 
Gi. These wave functions of the matter waves 
emitted from all slits of the gratings Go and Gi 
have been written out in Eqs. (25) and (26). 



Note only, that the wave function |^i) contains 
extra parameters, integer Ki and real Vji. 
Figs. [l4|(a), 14(b), and [l4|(c) show interfer- 



ence patterns emergent behind gratings contain- 
ing only A^o = 4 and A^i = 5 slits. Slits of 
the grating Gi are simulated by the set of the 
Gaussian functions ([38]) with the parameter r]i 
equal to 1.5 and for different Ki\ (a) Ki = 1] 
(b) Ki = 4; and (c) Ki = 16. More rapid diver- 
gence of the rays near the slits is seen to arise in 
the case of large Ki . These divergent rays form 
finely ruled interference fringes arising prior to 
the first Talbot length z^, see Fig. [T4|(c). 

Observe, that behind double the Talbot 
length, 2z^, the interference fringes become 



equivalent independently of choosing Ki. Quali- 
tative difference of intensities seen in Fig. [T4|(a), 
on the one hand, and in Figs. [T4|(b)-[T4|(c), on 
the other hand, is due to presence of the factor 
{2/tiY^'^ /rji ^ 0.5. This factor provides a good 



coincidence of a top plateau of the function (38) 
with that of the step function for Xi > 1. 
Whereas in case oi Ki = 1, height of the Gaus- 
sian function becomes twice as little. In partic- 
ular, the intensity in Fig. fl4|(a) is also twice as 



little than that in Fig. |4](a) for the same reason. 

One can observe a surprising phenomenon 
nearby the slits of the grating Gi as Ki in- 
creases - the larger Ki the better. The beam 
outgoing from a slit, first, goes through focus- 
ing. After it passes a smallest diameter (it over- 
comes so called the beam waist [32j) the beam 
begins to diverge. A detail pattern of this phe- 
nomenon in the vicinity of the central slit of Gi , 
see Fig. [l4|(c) , is demonstrated in Fig. [l5](a) . The 
outgoing intensity pattern shown in a range from 
O.bZrj. to O.GZrj. is seen to have a complex orga- 
nization. A tongue-like "jet" outgoing from the 
slit, first, undergoes squeezing. Dark patches are 
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Ki=16 



(a) 



3/2-Zj 



FIG. 14: Density distribution pattern p{x,z) in the near- 
field region z G (0, 1.5zt) = (0,0.15) m, No = 4, Ni = 5, de 
Broglie wavelength Ads = 5 pm, Talbot length zt = 0.1 m. 
Distance from Go to the source is Zs = —0.5 m, Xs = /xm : 
(a) Ki - 1; (b) Ki = 4; (c) Ki = 16. For ah cases ryi = 1.5. 



places where the intensity reaches maximal val- 
ues. The most maximal value is reached where 
the "jet" goes through the beam waist. Well or- 




0.52Zt OMZt 0.56Zt 

^1=64 



0.58Zt 0.6Zt 

(b) 



(a) (b) 



0.52Zj 



FIG. 15: Cutting out radiation from the central slit of Gi 
shown in Fig. 14 'c): (a) reproduced for case of Ki = 16; 
(b) reproduced for case of Ki = 64. Blue patches near the 
slit are most higher intensities of the radiation. 




x,//m 0.1 0.2 

FIG. 16: Profiles of the probability density p{x, z) within 
cross-sections (a) z = 0.52;^; (b) z = 0.5132;^; (c) z = 0.552;^ 
of the interference pattern shown in Fig. |15| 



ganized small dark patches seen in the figure pre- 
cede the most dark central patch. These small 
dark patches originate from center of the slit. 
And as they approach towards the most dark 
central patch, weak divergent rays can be seen 
are radiated far aside. More detailed jet outgo- 
ing pattern simulated for Ki = 64 is shown in 
Fig.|l5|(b). 

Fig. [16] shows some profiles of the probabil- 
ity density p{x^z) taken from the interference 
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pattern, that is obtained at given Ki = 64 and 



r]i = 1.5, Fig. 15(b). Integrals along the profiles 



give almost equal values. Insignificant discrepan- 
cies are conditioned by dissipation of the matter 
wave beyond an integration interval. 

The profiles (a), (b), and (c) depicted in 
Fig.[T6]are captured from cross-sections, that are 
marked by lines (a), (b), and (c) in Fig. [l5](b). 
They are at z = 0.5z^, z = 0.513z^, and z = 
0.55z^, respectively. The profile (a) abuts on 



the slit screen. We can see in Fig. [16] that the 
profile (b) clearly demonstrates against the back- 
ground of the profile (a) a squeezing of the beam 
together with increasing its intensity along the 
center. Next, the beam begins to disperse as 
distance from the slit increases, see, for exam- 
ple, the profile (c) in the same figure. 

C. Emergence of the focusing spot 

Emergence of the focusing spot with increas- 



ing number of the Gaussian functions in (38), 
i.e., with increasing Ki is an astonishing phe- 
nomenon at shaping slits with more hard edges. 
An effect of squeezing beam right after the slit 
can be understood at choosing of the slits that 
are approximated by the curve ([38]) with small 
Ki and r]i. Let us study such an effect (a) with 
increasing Ki at fixed r]i; and (b) with increasing 
rji at fixed Ki. 




FIG. 17: Difference Ap{x) of two profiles of the probability 
density p(x, z) given on the slit {za = O.Sz^) and in the cross- 
section of the beam waist (zb = 0.513^.^) as a function of Ki. 
The real parameter 771 is equal to 1.5. 

Let US consider, first, emergence of the focus- 



ing spot with increasing integer Ki from 1 to 
16 at fixed rji = 1.5. For this aim we pro- 
pose to evaluate a difference of two profiles of 
the probability density p{x^ z) given on the slit 
{za = O.bZrj.) and in the cross-section of the beam 



waist (zh = 0.513z^ 



Ap{x) = p{x, Zb) - p{x, Za). 



(52) 



As seen in Fig. [16] this difference can have a large 
positive hill in center of the beam bounded by 
wells from the both sides. All lies within a space 
of the slit. Fig. [17] shows a series of the dif- 



ferences /\p{x) to be calculated for Ki ranging 
from 1 to 16 and at fixed vji = 1.5. One can 
see, given Ki equal to 1 there is no Ap{x) dif- 
ferent from zero. As Ki increases from 1 to 16 
Ap{x) grows quickly enough and at Ki = 16 its 
peak reaches almost 0.1. It is the focusing spot 
bounded by the wells from the both sides. A 
depth of the wells is about —0.05 and they lie 
between edges of the slit, between —bi and 61. 



The difference (52) brings to light the focusing 



effect of outgoing from the slit the beam well 
enough. 

Let us now consider shaping the focusing spot 
with increasing rji but at fixed Ki. For defi- 
niteness, we will choose integer Ki = 7 and 
real rji ranging 0.2, 0.5, and 0.8, all are smaller 
than 1, see Fig. 18, Integrals of the approxi- 
mating curves are equal to a square of the step 
function for all cases. 

It is seen that the first curve. Fig. [T8](a), sim- 
ulates, in fact, seven slits situated within the in- 
terval [—6, +b]. Let 261 ^ 150 nm, then a period 
of such a finely ruled grating is about 20 nm. It 
corresponds to about 50 atoms situated between 
finely cut slits. The Talbot length, z^ = 0.1 m, 
earlier computed represents a very large length 



as against a Talbot length, z 



10- 



m. 



for 



the finely ruled grating having short intervals 
between the finely cut slits. We may observe 
in that connection series of principal diffraction 
maxima from this finely cut grating. These max- 
ima are partitioned by {Ki — 2) = 5 subsidiary 
maxima [29j . All these rays are divergent in the 
far-field region of the finely ruled grating which, 
in turn, is situated within the near-field region 
of the grating Gi . 
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FIG. 18: Approximation of the step function by set of the 
Gaussian functions presented in ( 38 ) with integer K — 7 and 
real (a) rj = 0.2; (b) ry = 0.5; (c) ry = 0.8. 



The above described picture is shown in 
Fig. [T9|(a). One can see, that a central princi- 
pal ray, diverging forward from miniature Talbot 
clasters packed by a triangle-like manner, has a 
narrow width before it will begin to diverge fur- 
ther. The width is positioned within z ^ 0.512z^ 
to z ^ 0.524z^. Positions of these two points in 



Fig. [19] are marked by vertical lines. 

Let us look on the second and the third curves 
in Figs. 18 [^b) and[T8](c), that have been simu- 
lated at Tyi = 0.5 and r]i = 0.8, respectively. We 
see, first, because of emergence of a pedestal un- 



der the function ( 38 ) at increasing the parameter 
77, intensity of a radiation from the slits grows 
up considerably. With increasing the parame- 
ter T] the function has a trend to approximate 
to the step function. Observe that, intensity of 
the radiation from the slit becomes the larger 



the more precise the function (38) approximates 
the step function. At increasing the parame- 
ter T] the triangle-like pattern near the slit is 








FIG. 19: Density distribution pattern p{x, z) right after the 
grating Gi\ z ^ (0.5^^^,0.54^^) and x G ( — 125,125) nm: 
(a) T] = 0.2; (b) T] = 0.5; (c) ry = 0.8; (d) ry = 1.1. i^i = 7 and 
A^o — 2^ Ni — 1. De Broghe wavelength Ads = 5 pm, width 
of the sht 26 1 = 150 nm, and Talbot length z^ = 0.1 m. 



seen to dissolve in the main ray. Together with 
that, the focusing spot is formed after the point 
z ^ 0.512z^. As the approximating function 
approximates to the step function, the focusing 
spot of the slit's beam becomes clear apparent 
within the interval from 0.512z^ to 0.524z^, see 



Figs. 19(b) and[l9Kc). This is an interval where 
rays formative the central principal maximum 
undergo squeezing in front of subsequent diver- 



gence [29j. Fig. 19(d) shows a formed beam right 
away after the slit in case of the parameter 77 > 1. 
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VI. CONCLUSION 

By utilizing the path integral method, we 
have computed the interference pattern from two 
gratings, placed in consecutive order along a par- 
ticle beam. A wave function describing the in- 
terference is found by summing all possible tra- 
jectories of the particles passing through slits in 
the first grating and then in the second grating. 
This powerful method permits to describe the in- 
terference eflFects of the matter waves in details. 

It is instructive to attract attention to manifes- 
tation of wave-particle duality nature at describ- 
ing the wave processes by this method. At first, 
we begin to consider particle paths going from a 
source to a detector through all possible interme- 
diate points. Ensemble of all possible paths gives 
description of effect of propagation of a matter 
wave. Wave fronts are represented by equiphase 
surfaces secant the trajectory bundles. Observe, 
that these waves are those that underlie in the 
Huygens-Fresnel principle [37] . 

The wave function found by the path integral 
method gives clear picture of the interference 
both between the gratings and behind the second 
grating. The most impressive observation of the 
interference effect is emergence of the Talbot- 
like carpets when illuminating the gratings by 
distributed coherent sources that are situated 
near the gratings, see Fig. [5](a). The carpets 
are smeared out at illumination of the gratings 



by incoherent matter waves, Figs.[5](b) and[5](c). 
Emergence of pedestals, supporting the inter- 
ference fringes, and decreasing the visibility are 
conditioned by breaking coherence of the beam 
that results in smearing of the interference pat- 
tern. It confirms results obtained at experimen- 
tal observation of the interference fringes fulfilled 
on molecular beams [101 1121 SHI • 

The grating, prepared with more hard-edged 
slits, as was shown corrects particle fiows in the 
vicinity of the slits the more powerfully than the 
grating with the fuzzy edged slits. An aston- 
ishing effect of the more hard-edged slits is that 
they shape focusing spots of the particle beam 
just after the slits, see the density distributions 
in Figs. [TSl and [TSl Next, they disperse along the 
beam with producing ripples. The situation is 
the same as rays initially converging within a fo- 
cus and diverging after. The focusing spot looks 
as a tongue-like jet passing through the beam 
waist. In particular, Nye [39j has shown that an 
electromagnetic monochromatic plane wave, in- 
cident on a perfectly conducting screen of vanish- 
ingly small thickness that contains an infinitely 
long slit, reproduces behind this slit an analo- 
gous tongue-like EM jet. Increasing EM ampli- 
tude right after the slit is shown in this case as 
well. One may suppose, that the phenomenon of 
focusing by the slits having more hard edges is 
typical for many wave processes. 
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